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Abstract 

We consider 0(iV)-symmetric potentials with a logarithmic singularity in the second field deriva- 
tive. This class includes BCS and Gross Neveu potentials. Formally, the exact renormalization 
group equation for the Legendre transform of these potentials seems to have ill-defined initial con- 
ditions. We show that the renormalization group equation for the local potential has well-defined 
initial conditions and that the logarithmic singularity is smoothed rapidly in the flow. Our analysis 
also provides an efficient method for numerical studies. 
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I. INTRODUCTION 



In quantum field theory and quantum statistical mechanics, bosonic O(N) models orig- 
inate naturally from microscopic fermionic models as the effective low-energy models for 
order parameter fields, like Cooper pairs or spin operators. As such, they play a central 
role in the analysis of symmetry-breaking phenomena. Technically, they arise via the in- 
troduction of auxiliary boson fields <fi coupling to composite fermion fields with a Gaussian 
integral (Hubbard-Stratonovich transformation). The correlation function of the composite 
fermionic "order parameter fields" can then be expressed as functions of the correlations 
of the 0, and the integration over the fermionic fields yields an new action Go(0) f° r the 
bosonic fields 4>. There are situations where the resulting action Go is not localized enough, 
so that the fermionic degrees of freedom need to be kept even at the lowest scales, but there 
is a large class of models where studying G is justified at low enough energies. 

On the mean-field level, a nonvanishing expectation value of the Hubbard-Stratonovich 
field <p signals symmetry breaking. In the full theory, fluctuations need to be taken into 
account (and can strongly change or even invalidate the mean-field result). 

The functional renormalization group (RG)^^ 1 ^ is a very useful tool for studying such 
fluctuation effects: it defines a flow of effective actions G s , with initial condition given by 
the potential Go, as a function of a scale parameter s > 0. In the typical application, s is 
related to some energy, length or temperature scaled that labels which degrees of freedom 
are incorporated. Here we have taken the convention that the energy scale is a decreasing 
function of s (or the length scale is increasing in s). We remark in passing that the RG 
method is flexible enough to allow for widely varying choices of s. 

More and more fluctuation effects are incorporated as s increases, and the full generating 
function for the correlations is obtained for s — > oo. The existence of this limit is not 
obvious. Indeed, control over this limit can be considered the solution of the model, i.e. the 
construction of a particular model of quantum field theory or statistical mechanics. 

There are several different implementations of the RG ideap"^, all of which are equivalent 
on a general level, but each with their proper merits and drawbacks when doing analysis 
and making approximations. The RG differential equation for the generating functional 
T s for the one-particle irreducible (1PI) vertices of O(N) models has had success in a wide 
range of applications, see Ref. O for a review. 
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It is important to note that the often-studied case of smooth initial potentials, e.g. 4 
potentials, does not really correspond to a model derived from integrating out the fermionic 
degrees of freedom. In many examples, the second field derivative of the boson potential 
contains a logarithmic singularity for small fields. The most prominent example is the BCS 
theory of superconductivity^, where the order parameter describes the superconducting 
gap. The same logarithm in the second field derivative can be seen in the Gross Neveu 
model 1 ^ and is relevant in the study of mass generation and chiral symmetry breaking in the 
two-dimensional situation, where the model is perturbatively ultraviolet renormalizable. 

This singularity in the effective potential cannot be regarded as a physically irrelevant 
detail because it implies the persistence of a symmetry-broken solution down to arbitrarily 
small values of the interaction strength. Indeed, all the familiar formulas of BCS theory 
would change if the potential were nonsingular. (Other features of the fermionic effective 
potential are not well-described by a 4 type potential either, as discussed below.) 

In this paper, we discuss the role of such initial singularities in the RG flow. Let W 
be the generating functional of the connected correlation functions. Following Ref. Q we 
set up the RG flow by multiplying the integrand of the functional integral for W with a 
regularizing Gaussian exponential with covariance c s = R~ x to obtain a scale dependent 
generating functional W s . Here s is the RG scale, which runs from zero to infinity, and R s 
is a regulator function chosen such that in the limit s — > 0, R s — > oo, so that all fluctuations 
are suppressed at the beginning, and the generating function for the amputated correlation 
functions is equal to the initial action Go- In the opposite limit s — > oo, R s — > 0, so 
that the regulator disappears and formally, the full generating function for the correlations 
is recovered (as mentioned above, it is nontrivial to show that this limit really exists). 
Taking the Legendre transform of the logarithm of the partition function, subtracting the 
regulating Gaussian exponent, and differentiating, we obtain the 1PI flow equation of a 
modified Legendre transform 2 - 



r.M = ^Tr 



In comparison with Ref. 3, Eq. ([T]) originates from a normalized partition function, that 
is, a term ^TrRgR' 1 is subtracted here. This functional equation is exact, but in most 
physically interesting models, the functional t— > T s [0] has to be approximated for a direct 
computation. 
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There are two common approximations for the functional r s [</>]. First, T s can be expanded 
in powers of the fields (ft an d truncated at some finite even order. If the local potential of 
T s contains logarithmic terms in (ft, an expansion around (ft = is obviously not possible. 
However, as discussed, the logarithm for small fields in the second field derivative generically 
ensures a nonvanishing mean field solution M . F .. By changing the expansion point to 4> m .f., 
one can avoid the logarithm in an expansion in (ft — M F . However, this expansion can then 
converge at most for \(ft\ < |</> m .f.|, which is very small for weak interactions. Even worse, for 
the BCS-model we find that the coefficient of {(ft — </> m .f.) 4 has a negative eigenvalue in the 
radial mode. Therefore, requiring stability of the functional integral, a 4 truncation is not 
feasible in this case. This problem is not cured by including the six-point function or by a 
naive separation of small and large fields. It is, of course, merely a problem of the expansion 
in powers of (ft — (ft M , F ., since the potential is bounded below. - We note in passing that 
potentials obtained by the Hubbard-Stratonovich transformation and fermionic integration 
also do not grow like \(ft\ at large \(ft\ but rather like \(ft\ 2 since the logarithm of the fermionic 
determinant grows only linearly in \(ft\ at large \(ft\. 

The other often-used approximation is a derivative or gradient expansion^^^i. While 
it is not yet clear under which circumstances such expansions are asymptotic^, they have 
been applied successfully to a variety of physical problems in a renormalization context, see 
Ref. I2I and the references therein. A naive application of the derivative expansion meets an 
ultraviolet problem for the case of fields (ft originating from a Hubbard-Stratonovitch trans- 
formation, because the fermion loops determining the action Go vanish at large Matsubara 
frequencies. Therefore, a time derivative term is never really there to smoothen the short- 
time fluctuations, i.e. the propagator for (ft has no decay at large frequencies. Ultraviolet 
divergences are only prevented by the decay of the higher vertices of the initial action in 
these frequencies, i.e. the decay of the vertices generated by the ^-dependent terms in ^fi 5 -. 
When the initial action is the result of an integration where the high-frequency modes are 
integrated over, e.g. in a fermionic representation, this ultraviolet problem is absent. The 
fermionic integration over high-frequency modes can be done by convergent perturbation 
theory.— 

A further problem is that the status of ([1]) becomes unclear in the limit s — > if the 
second field derivative of the initial interaction potential contains a singular term, such as 
^T#[0 C ] ~ ln0^, when evaluated at a constant field (ft c . Certainly, if one tried to replace 
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by in the inverse in Eq. ([T]), one would end up with a singularity at some small, 
c -dependent s. 

In the present note we show that this problem is not really there, due to the smoothing 
properties of the RG flow, which become evident when regarding the flow of the connected, 
amputated functions instead of the 1PI vertex functions. We show that the generating 
function for the connected functions is smooth at any s > and use this to give estimates 
on the Legendre transform that imply smoothness of T s in (ft for any positive s. We apply 
this in two ways. First, we can overcome the problem of the seemingly ill-defined initial 
condition simply by the semigroup property of the RG: performing the fluctuation integral 
with covariance c e as a Gaussian convolution for W e and respectively G e , and then Legendre- 
transforming, gives a new, smooth, initial condition T e for the generating function of the 
1PI vertices. It turns out that ~ logc £ , so that c £ ^f vanishes as e — > 0, and hence 
there is no singularity in the inverse in Eq. (CQ). Second, we use these estimates to show that 
the differential equation for the 1PI vertices holds for any s > 0, and we give the asymptotic 
behaviour of the solution for small s > 0. As one would expect, the deviation from the 
initial condition To is nonuniform in (ft, which explains the absence of the above-mentioned 
singularity: at any s > 0, one can choose (ft so small that To[(ft] is not a good approximation 



Thus the physically important logarithmic singularities in the initial condition for the 
potential do not present any conceptual problem for the functional RG, and our method 
also provides a practical method to treat such initial conditions, also in the 1PI scheme. For 
simplicity of presentation, we concentrate here on reduced O(N) models, that is, only on 
the local potential. The field theoretical methods and the estimates we use allow, however, 
generalize to the full model: the smoothing property of the Gaussian convolution also holds 
for infinite-dimensional Gaussian integrals, and the strong decay properties imposed by the 
RG regulator function at the beginning of the flow justify perturbation theory. In particular, 
the generalization to include the second order of a derivative expansion is straightforward. 

Let (ft = (0i, . . . (ftjsf) G M. N be a constant field, that is, a vector with N components. For 
H E R N let ((ft,H) = YlLxfa 11 " and denote ^ = (<f>,4>)- We consider a reduced O(N) 
model with the generating function for the connected correlations 



for r s [4 



W s (G Q ,H)=\n 



1 




^exp [ -^--G ((f>) + ((f>,H) 



(2) 



(2tts) n /' 
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The external field H couples linearly to 0. The scale dependence c s = sIn with s G [0, oo) is 
already included in the definition. This particular choice of scale dependence is not essential 
for the calculations; it is chosen for convenience only. The potential Gq is 0(./V")-symmetric, 
so that it can be written as Gq(4>) = Vo(p) with p = \<fi 2 - We assume that Vq is smooth 
away from p = and that for large p, Vg'(p) > const. > 0. For small p, we assume 

V (p) = Vb(0) + v lP \np + v 2 p + K{p). (3) 

Here V\ > and the remainder term 1Z satisfies 72.(0) = TV (0) = 0, and there is a constant 
K > such that \Tt"(p)\ < K p~ a with a < 1. With these assumptions, the function 
exp(— G o (0) + {H,4>)) is integrable uniformly in if, hence the limit s — * oo of ([2]) exists by 
the dominated convergence theorem. 

An important example satisfying these hypotheses is the mean-field potential of the BCS 
model. This is the case N = 2 and 

V (p) = 9 -- JdE u(E)y/E* + p, (4) 

if the density of states u(E) is regular at the Fermi level E = 0. Here — g is the coupling 
constant in front of the Cooper pair interaction term. The logarithm in (jHJ) is really there, 
i.e. vi > 0, if i/(0) ^ 0. 

For notational simplicity we have used a unit volume here. In general, the exponent is 
given by QVq, where Q denotes the volume, which is taken to infinity in the thermodynamic 
limit. In this limit, Eq. (j4j) becomes exact for the reduced BCS modet^. In presence of 
Q, the factor |0| 7V ™ 1 in the integration measure, d N <p ~ |0| Ar_1 d|0| d^^" 1 ^;, where d^" 1 ^; is 
the integration measure of the (N — 1) dimensional sphere, is not relevant for the following 
discussion, because all other parts of the exponent get multiplied by Q. 

The effective potential T <.(</>) = j s ((j)) — f-, where 7 S is the Legendre transform of W s , 
is again 0(AQ-symmetric and we write T s (0) = U s (p) (recall that p = \<fi 2 )- Denoting 
differentiation with respect to the scale s by a dot and differentiation with respect to p by 
a prime we obtain the RG equation 
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(N-1)U' S , U> s + 2pU» 



(5) 



l + sU> 1 + s [U' s + 2pU'j] 

for the effective (local) potential^, which can also be derived by inserting constant fields 
in Eq. (pQ). In this sense U s is the lowest order of a derivative expansion. Formally, the 



initial condition is posed in the limit s — > 0, where Uo(p) = Vo(p), which seems to lead 
to the vanishing-denominator-problem discussed before because U (p) = Vihxp + V\ + v-i + 
O^P 1 °) ~> ~°° as p — > 0. Of course, W s is convex by Jensen's inequality, and hence the 
Legendre transform cannot diverge at any finite p. In the following we show the more specific 
statement that, due to the smoothing effects of the RG transformation, the denominators 
are strictly positive, and we give sharp bounds for their behaviour as s — > 0. 

II. THE RG AS A SMOOTHING OPERATOR 

We introduce the effective action 



such that W s (Go, H) = | — <j s (£) with £ = sH. By 0(N) symmetry we can write G s (£) = 
V.(() with C = £ 2 / 2 - The structure of ® is 



where * denotes convolution and p s is the Gaussian measure with covariance s (the integral 
exists by the above-mentioned properties of Go). For s — > 0, p s tends to a Dirac measure, 
so the convolution gives e _Go in that limit. The convolution with a Gaussian measure is 
a standard example of a smoothing operator—, so this already implies that in spite of the 
singularities in derivatives of Gq, p s * e~ G ° is smooth, even analytic in for any s > 0. This 
can be seen explicitly from (p s */)(£) = / f(x)dp s (x — £), and understood in a physical 
analogy by noting that the RG flow defined in ([6l) is a heat flow with time parameter s, 
whose solution is smooth for any positive time s > 0. 

Therefore we can avoid the singular initial condition altogether by using the semigroup 
property^ of Gaussian integration: let e > 0, then for all s > e 





G.(0 = -ln {p s *e- Go )(0 



(7) 



G s {0 = - In (p s -e * e 



)(£)■ 



(8) 



Or in terms of the unamputated connected functions with a shifted scale 



W S (G E ,H) = W S+£ (G , 



s + e 



s 



H) + 



H 2 se 
2 s + e' 



(9) 
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for all s > 0, i.e. 

W S (G £ , H) — In J j^wf, exp [ - ^ - G e {<j>) + (0, )J • (10) 

We find lim™(iy s (G £ ,#) - W S (G 0} H)) = e^-, that is, the functions W S (G 0: H) and 
W s ((j e , H) coincide in the limit s — > oo up to an explicit term. The RG flow of the (modified) 
Legendre transform remains unchanged but the advantage is now that the initial condition 
of Eq. is given by G e (<j)) = V e (p), which is smooth. In the remainder of this section we 
compute V £ and give bounds on its derivatives. V e has no logarithmic divergence in for 
arbitrarily small e > 0, and it provides a well-defined starting point for integrating (jHJ). 

To begin, we collect some properties of Vq that follow from ([3]) and the assumptions on 
the remainder term TZ stated there, namely that, loosely speaking, the behaviour of Vq is 
that of V\p\np for small p. By our assumptions and integration in p, 

\tz"(p)\<^, |^(p)|<-^V-\ \n( P )\<J^ P ^ (ii) 

p a 1 — a 1 — a 

with a < 1. It follows immediately that 

\VM - v x lnp| < Vl + \v 2 \ + -^p 1 "*, (12) 

1 — a 

which is much smaller than |vilnp| for small enough p, and 

WW - v j\<^ (13) 

which is again much smaller than — for small enough p because a < 1. The properties 
of p h-» t>i In p and an easy approximation argument then imply that there is an interval 
(0, 2po] on which the derivative Vq of the initial potential is negative, the map p — > \Vq(p)\ 
is decreasing and the maps p — > p\Vq(p)\ , k = 1,2 are increasing. Moreover, on this 
interval |V^"(p)| < — where c" is a constant. In particular we can choose e so small that 
e\ Vo(e)| < 0.1. For reasons of brevity, we do not give the detailed values of the constants as 
functions of Vi,V2,K and a here. 

We split the analysis of G £ in two cases distinguished by the value of ( = y. 
Case 1: ( < £ < 1. We change integration variables to <f> = (0 + ^)/y/e, subtract V (0) 
in the exponential, and expand the exponential of Vo(0) — Vo(pe), where p = Then 
perturbation theory for small e yields 

1 



e -G«(f) _ e -Vo(0) 



l-e{v x \ne + v 2 )^(N + J 2 ) 



+ ev 1 T(J) + 0((e\ne) 2 ) , (14) 
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where J = R-sJe = G [0, y/2] and 



T(J) 



e 2 



IN, 



pin pe~ 



(15) 



(2tt)^/ 2 

The function T(J) and all its derivatives with respect to J are bounded on the interval 
J G [0, \/2\. Likewise the higher order terms and their derivatives with respect to J can be 
estimated. That is, although e can be arbitrarily small, V £ contains no logarithms of the 
field anymore. Additionally we obtain for the derivatives 



V^C)=vi\ne + v 2 + v 1 f 
V:'(0 = jf"(l) + O((\ne) 



(16) 



Case 2: e < ( < pq. We perform the integral (jBJ) by the saddle point method (because we 
are analyzing G s for s = e, s is substituted by e in (jHJ)). The stationarity condition for the 
negative exponent S(<f>) = + Vo ( \ ) in the integrand of Eq. (Ej) is 



OS 



f Vo 
e L 







(17) 



for all i. We first assume that there is a stationary point 0* and denote p* = (0*) 2 /2. Then 
( TTTjl implies 



2 = 2e 2 p*V \p*" 2 



(18) 



and 



P*U+^oV)) 2 



C- 



(19) 



The left hand side of ( jT9l) is monotonically increasing in p* G [e, 2po] by our hypotheses on 
the potential V . Thus a unique solution p* G [(, 2£] of (IT91 exists. There is no solution 
in the interval [0,e] since C > e. For larger fields there is no solution since V Q \p) becomes 
positive eventually, so that Eq. (JT9"j) would imply p* < ( < p , and because e is small. Given 
p*, the unique solution of (fT7|) is, by O(iV) invariance of Vo, 0* = 
single stationary point. By ( fT9l) . and because V^(p*) < 0, 



2p* |||. Thus S has a 



0<p-(<p* {2e\V '(p*)\ + eX(p*f) 



(20) 
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Because \Vq \ is decreasing and p* > ( > e, this implies 

< p* -C < P*V(2 + V) < %P*V (21) 



with r] = e|Vo(e:)|, hence 



<<p-<rSr,- < 22 > 



We thus have the estimate 



IWO-K'(C)I<(p*-C) sup IK'COI 

re[C,p*] 



c" c" 



<(p*-0-<(p*-0 

C p*(l-3rj) 

< C "J^_. (23) 
-l-3?7 



These bounds imply that all eigenvalues of the Hessian 



H *i = l^T = -M l + eV o(p)) + ^X(P) (24) 



"3 £ 



are positive and of order e~ x at 0*. Thus <fi* is the unique minimum of S and a standard 
saddle point analysis^ applies: all contributions from <fi not in a neighbourhood of the 
minimum are suppressed exponentially for small e, as are the corrections to the Gaussian 
integral around the saddle point. The Gaussian integral around the saddle point gives 
(2n) N l 2 D~ l l 2 , where D = det H. It gives only subleading contributions since the factor e N ^ 2 
in D~ x l 2 is canceled by the normalization factor s~ N / 2 = e~ N / 2 of Eq. (jH]). Therefore in case 
2, 

Vs(0 = V (0 + O(e) 
V:(0 = V^0 + O(elne) 
CK"(0 = Oil) . (25) 

Combining Eqs. ffl6l) and (1251) from the two cases, we find that in V £ (Q the logarithm of the 
field ( is replaced by the logarithm of maxje, (}. Therefore, the RG flow starting at s = e, 
and with with initial condition T e = V e , is well defined. 

As we have just shown, perturbation theory for small ( < e allows us to calculate V £ to 
arbitrary precision. This result can easily be extended to non-reduced models because in 
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general, the regularization c e provides an infrared regularization, which justifes perturbation 
theory for small enough e. For reasons of brevity, we have only outlined the saddle point 
argument that estimates the difference of V £ and Vq for ( > e. This argument can easily be 
made into a proof, and it also extends to the non-reduced situation, again by noting that 
the infrared regularization together with the smallness of e provide rigorous control over the 
saddle point expansion. 

III. THE RG DIFFERENTIAL EQUATION AT SMALL s AND <j> 

Shifting the initial condition of the flow as described in the last section is an exact 
procedure and approximations become necessary only for the calculation of the new initial 
condition (at least for non-reduced models). But the question remains whether one can 
find a less indirect way of showing that the RG equation (jSJ) for the local potential U s is 
well defined at all s > if the initial potential contains logarithmic terms. In this section 
we study the asymptotic solution of the RG equation for small RG scales s and small field 
squares p = \4> 2 . As explained below, the argument is not solely based on (jSJ), but requires 
the bounds derived in the last section as an a priori input. 

In a first step, we assume that the denominators and also pU'J in Eq. (jSJ) do not contribute 
to the leading asymptotic solution. Then the flow equation becomes a partial wave equation 
U s (p) = yU' s (p), which is solved by the backward propagating wave 

U s (p) = U (p+fs). (26) 

If we knew that Eq. (1261) also provides the asymptotic behaviour for the derivatives with 
respect to s and p, we could easily justify the assumptions we just made: the denominators 
for small p and s contribute only to order O(slns), and pU" is bounded by a constant for 
small p. However, asymptotic expressions cannot simply be differentiated, hence regularity 
of the derivatives of the local potential cannot be assured by this argument. The natural 
procedure starting from the RG equation would now be to differentiate Eq. ([5]) with respect 
to p. This allows to determine the asymptotic solution and to verify the above assumption 
for U' s , provided that a regularity assumption is made on U'J . Another differentiation allows 
to do the same for U", given a suitable hypothesis on U' s ", and so on. To avoid an infinite 
proliferation, it suffices to have a priori bounds for U' s (p) and pU'J(p) for small s and p. 
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We have already derived such bounds directly from the functional integral in the previous 
section, and use them now to prove the asymptotic correctness of ff26l) . 

For the Legendre transformation of W s we denote the inverse of the maps ^^(H) \—> <fii by 
Hi(4>) = %^ + y- Using W S (H) = | — V S (Q we find the connection between the derivatives 
of the effective action V s (() and the local potential U s (p) 

u M= ™ (27) 

where ((p) = s 2 H ^ is determined by £(p) = p/[l — sVg(C(p))] 2 ■ Combining the estimates 
of Vg(Q obtained in Eqs. (TTB1) and (]25l) we arrive at the estimate |V S '(C)| < cln(max{s, (}) 
for small s and ( and a constant c 6 R. Using Eq. ( 1271) this gives |t^(p) | < cms with another 
constant c. Similarly, |U,"(C)| < c(max{C, s}) -1 implies < cp _1 asymptotically for 

small p and s. Therefore, Eq. (|26|) is the asymptotic solution of the RG equation. 

IV. CONCLUSION 

We have shown that it is possible to apply the functional RG to initial conditions given 
by potentials with a logarithmic singularity in their second field derivative, because the RG 
flow smoothes out these logarithms sufficiently fast. One might think that a rapid change 
of the effective local potential near the singularity might cause numerical difficulites, but 
our arguments also provide a method to calculate the flow at small s efficiently and with 
arbitrary precision. 

We have restricted our analysis to reduced models to bring out the main points in a 
simple way, but it can be generalized to include the second order of a derivative expansion. 
For example, the Zq and Yq functions (see Ref. 2) for standard notation) diverge with p" 1 
and p~ 2 respectively for the BCS model. As shown here for the local potential, regularized 
functions Z e and Y e can be obtained by a derivative expansion of the effective action at scale 
e. Moreover, as explained above, the smoothing argument is completely general, that is, it 
can be used to prove a similar statement to the full theory. 

As already remarked in the beginning, potentials with singularities are not academic 
examples, but arise in important physical situations and have important effects. The results 
described here will therefore be useful in going beyond 4 -type approximations of these 
potentials, to obtain a more quantitative theory. 
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A natural question is whether our analysis also applies to more singular initial conditions. 
It is straightforward to extend our proofs to potentials Vq whose derivative diverges as a 
power of logp for p — > 0. This case includes, in particular, a (log/?) 2 singularity, which 
occurs in the study of superconductivity of two-dimensional Fermi systems with Van Hove 
singularities. 
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